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ABSTRACT! 

This thesis deals with the computer simulated, elastic link model 
of an automobile cam actuated valve train. Examined are the dynamic 
responses to such situations as: valve bounce, pushrod bounce, 
excessive cam speed, cam surface machining errors, and cam profiles. 
The model is adaptable to various types of mechanisms and lends itself 


well to the incorporation on a computer graphics display. 
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IT. INTRODUCTION 

The study of cam actuated valve train mechanisms has been undertaken 
in the last thirty years from many different approaches. The mechanism 
has often been dissected to observe various local responses with the 
remainder of the mechanism considered as external to the local point of 
interest. 

Same history of past work undertaken will help to demonstrate the 
usefulness of the results to be presented later in this thesis. 

Of basic interest to those dealing with valve train mechanisms 
is what role does the cam profile play on the over-all response of the 
system. Dudiey [1], [2] examined the design of cam profile by use of the 
valve output. This technique was carried out also by Thoren, Engemann, 
and Stoddart [3] who included experimental results based on observations 
Of various engines and cam designs. More recent analysis has been done 
by A. R. Johnson [4]. Johnson dealt with a mimerical applicaticn of 
Newton's Fundamental Interpolation Formula to look at a multi-degree of 
freedom system. Again the idea was to take the desired output and relate 
Hemet terms Of the input. The above references are significant for they 
reaffirm the importance of the output being a dominant factor to the design 
of a cam. 

Historically there has been a large amount of controversy over the 
best cam for design purposes. Three types of cams: parabolic, harmonic, 
and cycloidal often have been studied and compared to each cther in terms 


of displacement, velocity, acceleration, and jerk [5], [6], [7], [8], {9]. 


Polydyne cam profiles with their added flexibility have beceme: 
popular. They can easily be adapted to sets of conditions which may be 
imposed by the mechanism [1], [2], [4], [10], [11], [12], [13], [14]. 

The problem here is that the cam may be over designed, in other words, it 
may be too exactly defined to ke manufactureble. 

In addition to those mentioned above, references [15], [16], and [17] 
present other methods of obtaining a cam profile. As can be seen there is 
a wealth of information to the cam designer to aid him in obtaining a cam 
profile, provided he has the necessary input data. 

The ability to analyze a particular cam profile's effect on a system 
can be of great benefit. Experimental results are the best, but also the 
costliest. In addition, experimentation is best utilized if it can be 
compared with theoretical results. This can give indications of existing 
conditions which may not have been considered previously. 

The ability to produce a model of a mechanism provides the chance to 
obtain a balance of varied information for analysis. Models have the 
same deficiencies as the human brain, i.e., they consider only a portion 
of the existing conditions. They are often simplified and are progranmed 
to fit a situation as the person establishing the model sees it. This 
procedure can be less than perfect as can be pointed out by experimental 
results. On the other hand, it ee as a good guideline and usually can 
be handled, revised, and adjusted more easily than a laboratory set up. 
Models have been used extensively in the investigation of cam actuated 
systems. Most of the work has been done by the classical method of 


mechanics, i.e., a free body diagram with force analysis, [1], [2], [9], 


[18}, [19], [7], [20]. The technique used is to set up the free body diagran 






and solve the resulting differential equations. Graphical application 
to solve for acceleration, velocity and displacement polygons can also 
be used [21]. 

The above references usually employ a mechanism that behaves properly, 
i1.e., the cam and pushrod stay in contact and problems such as pushrod 
bounce, valve bounce, impact loads, cam machine surface errors are not 
considered, or are considered alone and not in connection with each other. 
Barkan and McGarrity [20] did consider impact loads on takeup of initial 
cam clearances, and R. C. Johnson [22] stucied the impact force in a 
system using *inite mass and infinite mass assumptions. First he calcu-~ 
lated the forces between two bodies of finite mass. Next he assumed one 
of the bodies, the driving element, to have an infinite mass and calculazed 
the force for this case. In this way he obtained a range of expected impact 
force values. Rothbart [7] in Chapter & and Neklutin {23}, [15] provide 
insight into the problem to be expected with vibrations internal toa 
mechanism. Mitchell (5), also considered the interaction of cam profiles 
with the resonance set up at the natural frequencies of the mechanisn. 
Rothbart, Chapter 10 [7] gives expected dynamic results of can surface 
machining error and the response to the system. 

One of the best and most comprehensive studies of valve train dynamic 
response is that of Barkan [19]. He established a valve train mechanism 
which he compared to experimental results. His medel was developed py 
the technique of lumped parameters. The lumped masses were connected by 
massless springs of rigidity determined by the particular ccmponents. 
Included in the mechanism was internal and external viscous Camping. Forces 
caused by internal combustion and by the valve spring were externally 


applied to the mechanism. 
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Differential equations used in solution of Barkan's mechanism were 
solved by an integration technique using successive time increment evalu- 
ations. This type of solution technique is also demonstrated in Rothbart 
[15] page 240, and Timoshenko and Young [24]. | 

R. C. Winfrey [25] adapted matrix algebra to the solution of the 
equation of motion developed in mechanism analysis. This, coupled with 
the calculation ability and speed of a digital computer, reduces the amount 
of labor involved in the arithmetic to a much more elementary level. The 
technique for obtaining the equations is that developed in Rubinstein's 
text [26]. 

Applying the tool, matrix algebra, to the model which will be described 
later gives the potential of studying many features of a mechanism response 
when subjected to various operating situations. The attempt is to produce 
a simplified, adaptable model which can exanune not only We ideal situations, 
but also the response to abnormal or undesirable circumstances. 

This thesis deals with situations which should not occur on normal 
valve train systems. This is done to impress the reader with the capa- 
bilities of the model. Examined will be impact loads, valve bounce, pushrod 
bounce, excessive cam speed, machining errors, and ene The model 
1s capable of giving a voluminous amount of data to pinpoint response 
analysis. It is adaptable to various types of mechanism analysis and lends 
itself well to the possibility of being incorporated on a computer graphics 


display. 


de: 














TI. ‘THEORETICAL MODEL AND APPLICATION 
OF EQUATIONS OF MOTION 


The model shown in Figure 1 is a simplified adaptation of the valve 
train system of an autamobile engine. It consists of a cam, pushrod or 
cam follower, rocker arm, valve stem, valve spring, and valve seat. The 
rocker arm pivots on the rocker. It is pinnec to the pushrod at one end 
and to the valve stem at the other end, thereby restricting all members 
to remain in contact. The valve spring is assumed to have a spring 
constant of 245 pounds per inch, is connected at the pin to the rocker 


arm and valve stem at one end and is fixed to ground at the other end. 


SPRING 
VALVE STEM 


VALVE SEAT (OPTIONAL) 
PUSHROD 





FIGURE 1. BASIC MODEL 
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The mechanism is acted on at one end by the cam in contact with the 
pushrod or cain follower. On the other end the valve may be acted on by 
the valve seat. The model is further simplified by each member having 
a constant circular cross section. 

The valve seat position is adjustable, dependent upon the desired 
output. It is identified by a dasned line in Figure 1 signifying that 
it is not always simulated in the system. For sav applications the 


valve is allowed to move freely, unimceded by tne valve seat. 


A. EQUATIONS OF MOTION 
The equations of motion,-: 

IMlig) + [Cliq}) + [Kl{q) = {F} (1) 
were adapted to the model. The coefficients of the mass, stiffness and 
‘damping matrices were determined by the finite element technique [26]. 

i, Syste Cooraimates 
The systen coordinates were established.as shown in Figure 2. 
It was assumed that these coordinates accounted for all relevant eneray. 
A linear model (1.e., small deformations) was also assumed. In this 
way tension end compression forces were all that were considered in the 
pushrod, valve stem, and spring, and only the bending forces in the 


rocker ann. 








FIGURE 2. SYSTEM COORDINATES 
ue Blement Coordinates 
Figure 3 gives the element breakdown of the system with the 
corresponding element coordinates. Again it is noted that elements 1, 
4 and 5 (1.e., the pushrod, spring, and valve stem) are treated as 
axial members and elements 2 and 3 (i.e., the rocker arm) are treated 
as bending members. The rocker arm was divided into two elements to 


tter represent the bending motion at the pivot. 





Yu, U, a, u, 
U eee ‘i iD 
: on t 
u. Uy us Oy fu, 2 


ne! 


FIGURE 3. ELEMENT COORDINATES 
3. System Mass and Stiffness Matrices 
The system coordinates, {gq}, are related to the element 
coordinates, tu} by a transforination matrix [8] as shown in Equation 2. 
{u} = [Bliq) 7 (2) 
[8] is formed by applying the condition of compatibility to the mechanism, 


and it is evaluated in Equation 3 for the model of Figures 2 and 3. 


ds 
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| | a 
Uv, Ye Oe a Oe eee q4 
U, = (0) idol emer: OM Pibeneeed (gaan 0) oe (0. (3) 
3 
u, O76 OE) OQ. (eels BO 
i 0 0 0 0 iL 0 0 
u 0 0 0 0 0 0 0 
1 A | 
+ RE eee 
U 6) @) 0 0 0 6) 1 
Sas 
: 0 0 0 1 0 O 0 
\ J 


As noted in Equation 3, [8] is partitioned into sub-matrices according 


to elements. Equation 3 may thus be rewritten in a more concise form, 


fu} y [8] 5 (eae 

fu}, [8} . 

{fu} , = [8] 5 (4) 
tu} , 8) , 

{fu}. [8] 


Utilizing the transformation matrix [8] and the procedure of Chapter 7 
of reference [26] the system mass and stiffness matrices, [M] and [K], 


were formed. 
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As shown in the reference, [8] relates the stiffness and mass properties 


of the elements to the stiffness and mass properties of the system. 


[K] (5) 


II 
fawn 

cc 
WD 

ies 
— 
~ 

is 
paemay 
WD 
es 


[M] (6) 


II 
i mo 
_—- 
WD 
is 
— 
S 
I 
py 
WD 
— 


[6] 4 is the submatrix of [8] corresponding to the ju element. Sim- 
larly, IK! ; and IMJ , are the i/4 clement stiffness and mass matrices. 


The element stiffness matrices are given by [25]. 








EB. A, 
KJ 45 aL iL (1 = 1,4,5) (7) 
1 
iL 1 
eid 
Be 3 
K].= 2 = 6L 6L 
=] sak L 3 a 
i: 
Aw, ie Zoi =GIr 
(3) 
6L -6L 4L,.2 Dee 
al sl 
{ 
lly, BGR « oc Jl sg (i = 2,3) 
a 1 1 BG 


and the element mass matrices are: 


M. 
IM}. = —+ 2 =i 
6 GeS Les) (9) 


1, 








And, 


M. 
el: | 
[M.] = Do 156 54 221, 
54 156 aeeet, 
Dor 13L. 4y,,2 
1 1 ai 
~13L. -22L oie 
3 1 





In the above equations, E. is the elastic modulus, A, is the cross 
Sectional area, L. is the length, and M, is the total mass of the “ae 
element. 

a. Stiffness Matrix 
Evaluating Equation 5 for [K] gives Equation ll (see page 19), 
b. Mass Matrix 
Matrix multiplication of Equation 6 yields for the systen 
mass matrix, Equation 12 (see page 20). 
4. Systen Damping Matrix 
The systen damping matrix 1S formed in two steps. First, the 
equation of motion without damping is uncoupled and is written in terms 
of modal coordinates. Next, by the use of the diagonal matrices of the 
normalized equation of motion, the damping matrix can be developed. 
| a. Equation of Motion in Terms of Modal Coordinates 
Rewriting the equation of motion, Equation 1, without 


damping yields: 


[MJ {fq} + [K]{q} = {F}. (13) 
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(12) 











Having the mass and stiffness matrices, the eigenvalue 


problem may be solved using Equation 14. 


(w2 [M] - [K]){q} = {0} (14) 
This gives a solution of seven eigenvalues, w (Gillen. 7) anG@ ese en 


corresponding eigenvectors or mode shapes, {¢}.. 


; Arranging the 


eigenvectors in colums, 


[o] = [{o}, {o} {e}, {o}.] (15) 


2 
gives the modal matrix, [¢], which may be used as a coordinate transfor- 
mation matrix to uncouple Equation 13. This gives the equation of motion 


in terms of the modal coordinates, 


IM {n} + [Kl {n} = (FY, (16) 


where {F} is the modal force vector, and [M] and [K] are the diagonal, 


or uncoupled, mass and stiffness matrices related to [M] and [K] by, 


tI 
c—_ 
«> 
— 
ie 
= 
ae 


[M] 
[K] = [¢]” [Kk] [¢] Jeo. (18) 


(17) 


{n} and {n} are a set of system modal coordinates related 


to the system coordinates {q) and {q} by, 


{q} = [¢]{n} (19) 
[o]{n} (20) 


I 


{q! 
b. Formation of [C] 


Having [M] and [K], fC) can be formed and [C] solved for hy 


Matrix manipulation and multiplication. 


elt 





—— 
<> —— Ss 
— = 





k, s 
s A: s 
The eigenvalues, Was equal ia with Kee and ms equal 


to the diagonal elements of [K]} and [M], respectively. The modal system 
is assumed to have viscous damping, so that the values on the diagonal 


Smeic) + Car are related to Kes and ms by, 


es 2c. Vag : (Sle, 2 <1) (21) 


where Ce is the damping ratio. The damping ratio, 7, of the fundamental 
mode was set equal to 0.5 and the other six modes of higher frequency 
were set equal to 0.1. 


The system damping matrix [C] can now be found from [€] as 


E@ulows : 


Thus, 
Ic) = [e7)~) te) fe)? (23) 
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III. ADAPTATION OF MODEL TO PROBLEM SOLUTION 

The model. established in Figure 1 and the equation of motion 
established in Section II can now be used to examine dynamic reactions 
or the model to applied situations. 

There are two coordinates of the model acted upon by external forces: 
coordinate one end coordinate seven. Coordinate one is the displacement 
Of the point on the pushrod which comes in contact with the cem, and 
coordinate seven is the dispiacamrent of the point on the valve which 
comes in contact with the valve seat. Depending on the velue of the force 
and displacerent at these two coordinates, the model can exist as one of: 
four possible configurations; (1) contact between the cam and pushrod ard 
no contact between the valve and valve seat; (2) no external contact at 
cam or valve seat; (3) contact between cam and pushrod and valve and valve 
Seat; and (4) contact between valve end valve seat and no contact between 
Saieaid pushrod. These four configurations will be referred to as con- 
figurations I through IV, and will be described in more detail shorcly. 

The model will fall into one of the above configurations at all times, 
but will remain in any configuracion only as long as the physical con- 
straints of the model are satisfied. Certain coordinate values, listed in 
Table I, are monitored to indicate when the model shifts from one confiogu~ 
ration to another. As an example, at coordinate one the force, Poy 
between the cam and pushrod must nave a negative value. The negative value 
is due to the choice of the direction of coordinate one and the fact that 
the cam can only push on the pushrod and cannot pull. Therefore, in the 


solution of F,, if the value changes frem a negative to a positive value 


neg 


ene mechanism is no longer acting within its comstraints and trieecme 
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TABLE. I 
CONFIGURATION RELATIONSHIP 


CONTACT 
CONFIGURATIONS COORDINATED KNOWNS 
I 1 Qj 19,19, = Cam 
(Pushrod) Profile 
Fl. = 0 i = Zao 
ie None FE. = QO 1= 1,2, 
nee Ley . Oy yy = Cam 
(Pushrod Profile 
and Valve) q5 = Valve Seat 
Position 
dq 1% = 0 
cs Ong epee 
a 
TV 7 q5 = Valve Seat 
(Valve) Position 


47197 = 9 


f. = 0 ae... - 


Al 


MONTTORED 
VALUES 


qyrFy 





mast shift to another configuration, in this case configuration II, 
which allows for no external forces. The value of Fy will remain equal 
to zero until the pushrod comes back in contact with the cam. 

The equations of motion for the four configurations differ slightly 
from one another. The rearrangement of the basic equation of motion, 
Equation 1, for each configuration is such as to take advantage of the 
known values of coordinate displacement, velocity, acceleration, and 
force for that particular configuration. The known values for each con- 
figuration are given in Table I. The adaption of the equation of moticn, 
Equation 1, to each configuration and the solution technique are given 
in Appendix A. ‘As can be seen in Appendix A, each configuration has a 
set of equations having the same form as Equation 1. This solution is 
done by the integration procedure described in Appendix B. An 


explanation of each configuration follows. 


A. CONFIGURATION I 

tars 1s the basic configuration. In this configuration, contact is 
between pushrod and cam only. Each of the examples studied later begins 
in this configuration. The force at coordinate one end the displacement 
at coordinate seven, are monitored to determine when the model is to 
Beane ato configuration II or IIi1.° | 

1. Force at Coordinate One 

Due to the selection of the direction of coordinate one, the 

force must be negative for any contact. A positive value for Foy desis 
cates a shift to configuration II, free of external forces. The value 


of the force at Fy is obtained from the first row of Equation l. 


Sa ag 


oe a ye a eo 


— 


Das: 


2. Displacement at Coordinate Seven 
The value of q-7 1s monitored in those problems in which a valve 
seat is considered. The valve seat is simlated as having a particular 
position and when interference between the valve and valve seat occurs 
the model shifts to configuration III, contact at both the cam and valve 


seat. 


B. CONFIGURATION II 
In this section the model is in the free or "floating" state. There 
are no external forces on the model and it moves freely based on its 
natural frequencies and its initial conditions when it assumed this con- 
figuration. The valve spring acts on the system as an internal member 
to bring the system back in contact with the valve seat or cam. The 
monitoring of Gy a i5 indicates wnich contact happens first. 
1. Displacement at Coordinate One 
qy is monitored against the position of the cam. When inter- 
ference occurs between the cam and pushrod, the model shifts into 
configuration I, contact between cam and pushrod. 
2. Displacement at Coordinate Seven 
d4 is monitored against the position of the valve seat and when 


it exceeds the valve seat position the model shifts to configuration IV, 


contact between valve and valve seat. 


C. CONFIGURATION III 
infeeconticguration Iii, both the Cam and valve seat are in contact. 
This configuration may be thought of as a cross-over configuration. The 


model will not remain in this configuration as long as it does in the 
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Gther three configurations. The model enters configuration III from 
configuration I or IV. In configuration I, the model has been in con- 
tact with the cam and in configuration IV, tne model has been in contact 
with the valve seat. 

Since the model is elastic, it will shortly lose contact at one of 
the two contact points. Viewing Figure 1] it can be seen that when the 
valve 1S in contact with its seat, and the cam comes in contact with the 
meenied, the end result will be to quickly push the valve away from the 
seat, and hence transfer to configuration I. When the pushred is 
rumning on the cam and the valve comes in contact with the seat, the end 
result will be to 11ft the pushrod away fron the cam and transfer to 
configuration IV. 


and F_ are monitored. 


The two external force values, Fy 4 


ieeeerce at Coordinate One 
The force at coordinate one iS monitored in the same manner as 
configuration I. When Fy shifts from a negative to a positive value, 
the pushrod loses contact with the cam and the model shifts to con 
configuration IV. 
2. Force at Coordinate Seven 


The force at coordinate seven is monitored by solving the last 


feoeon Equation l. 


7 . 
na og ate * 794 * H(i! ap 


Due to the selection of the direction of coordinate seven, when the value 
OL PF. changes from positive to negative, contact at coordinate seven has 


been lost and the model shifts to configuration I. 


Zi 





Dee CONFIGURATION IV 
Configuration IV indicates contact at the valve and valve seat only. 
The values Wy and Fo are monitored in the same manner as previously 
described and indicate when the model shifts to configuration III or II 
respectively. 
Figure 4 gives a schematic of the four configurations, their 
relationship to each other, and the coordinate values monitored between 


any two configurations. 
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CONFIGURATION I 
Contact 
Cam-Pushrod 


G9 
CONFIGURATION IIT ‘ CONFIGURATION IT 
Memtact Both No External 
Cam—-Pushrod Concace 


alve-Valve Seat 


CONFIGURATION IV 
Contact 
Valve-Valve Seat 


Values on diagonal lines are monitored to determine when to shift to 


another configuration. Arrow indicates direction of shift. 


FIGURE 4. INTERRELATIONSHTP OF CONFIGURATIONS 
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IV. IMPACT FORCE ANALYSIS 
A. BACKGROUND 

Impact in the model occurs when the cam and pushrod or valve and 
valve seat come together. Impact forces are difficult to predict accu- 
rately and most often are studied from an elemental standpoint, i.e., 
two idealized masses impacting with smooth surfaces and relatively 
Simplified geometry. 

Background investigation to determine the manner of handling impact 
forces in model analysis does not yield much useful information for this 
thesis. The cam was assumed to be an infinite mass with its velocity 
unchanged by impact. As was explained in Chapter I, this assumption was 
one of the two used by R. C. Johnson [22] to predict a range of values 
for impact. Timoshenko and Goodier [27] also studied colliding members. 
The above references dealt basically with two masses and did not consider 
an elastic mechanism with interacting members. Of prime interest was the 
local deformation and stresses produced at point of impact, whereas in 
this thesis, deformation of the entire system is taken into account 
rather than local deformation. 

Barkan and McGarrity [20] mention impact forces in a valve train 
system, but only from the standpoint of the takeup of initial clearance 
between cam and pushrod. Once their system began to operate there was 
no later separation. 

Dubowsky and Freudenstein, [28], [29], also deal with the problem of 
Clearances. in mechanical systems. They developed a dynamic model for 
their study and examined the system response. Their study deals with 
clearances at pins and points of connection and does not deal with tha 


effect of large impact forces on a system. 
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Thoren, Engemann, and Stoddart [3] did consider valve bounce as 
related to the design of cam contours. They conducted laboratory studies 
of different types of valve trains and with different cam profiles. They 
Gid not set up an analytical model for a valve train, and did not consider 
the effects of such situations as pushrod and cam separation and impact. 
In addition, they were mainly interested in valve displacements and did 


Noe measure contact forces. 


B. MODEL IMPACT FORCE CALCULATIONS 

The integration technique of Appendix B used in the solution of the 
equation of motion finds displacement, velocity, and acceleration over 
successive time ‘steps, AT. Figures 5, 6, and 7 give idealized curves of 
g and g for a particular coordinate which undergoes a sharp change in 
its displacement profile, gq. This produces a step input to the velocity 
curve and an impulse to the acceleration profile. Figures 8, 9, and 10 
are the modifications to the idealized case which are necessary to 
account for the modeling limitation of Appendix B. Changes can only 
Secu curing on integration interval, AT. 

Figure 8 is identical to Figure 5. Differentiation of the displace~ 
ment yields Figure 9 for the velocity curve. As shown it takes one time 
step to reach the step input value in Figure 6. Differentiation to 
obtain acceleration yields Figure 10. It takes two time intervals to 
account for the impulse in Figure 7 and reduce the acceleration curve 
again to zero. The value of acceleration at time t. 1s given by the 


2 
equation, 


g el Wes (28) 
aS 


ol 








time 





FIGURE 5. Ideal Displacement Curve 


time 


FIGURE 6. Ideal Velocity Curve 


- time 





FIGURE 7. Ideal Acceleration Curve 
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FIGURE 8. Model Displacement Curve 





FIGURE 9. Model Velocity Curve 


7 ENA 


Eiieo 3 
FIGURE 10. Model Acceleration Curve 
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Note also, that due to the value of ora, which will usualiy be 
2 


fairly high since AT is smali, the value of F the force at the 


t f 
Z 
coordinate in question, can be significantly higher than Fy or F 


This is to be expected at the instant of impact. 
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VR SUpIS 
Figures ll and 12 are the numerical values of the [M] and [K] matrices 
respectively. These were obtained by first substituting the values of 
Table II into Equations li and 12 to get the mass and stiffness matrices 


[M] and [Kk]. 


MASS MATRIX [M] 


15 -7.8 0.0 0.0 0.0 0.0 0.0 

-7.8 Deut ilgule 0.0 0.0 -.6 0.0 
0.0 ee: 2 00.0 Oeie = Leal 0.0 
0.0 00.0 0.0 20.1 -2.4 ed -2.6 | xl07-° 
0.0 0.0 0.0 -2.4 6 -.5 0.0 
0.0 -.6 - 1 tipi es 9 0.0 

q ae 0.0 0.0 Da ee On (0 5.2 | 
FIGURE 11 


STIFFNESS MATRIX [K] 


2.4 2.4 0.0 0.0 0.0 0.0 0.0 
me 25.2 11.5 © 0-0 90.0. 11,5 0.0 
mor 11.5 7.6 0.0 0.0 3.8 0.0 
0.0 0.0 0.0 WG Sell 3. Bol Fal 105 
0.0 0.0 0.0 Sol svi, “OG 0.0 
moe 11S) 3.8 Sc SS OG 0.0 
0.0 0.0 0.0 Jala On O08 fa 
FIGURE 12 





TABLE IT 


Blement Length {L) Inches Diameter (D) Inches 
uy O30 ss 
DB ea 6 
3 Ac 6 
4 Son 0 07 
5 a 3 


FOR ALL MEMBERS 
Elastic Modulus = 30 x 10® lb/inch? 


Specific Weight = .283. 1b/inch? 
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With the mass and stiffness matrices formed, the eigenvalue problem 


was solved and the damping matrix [C] obtained, Figure 13. 


DAMPING MATRIX [C] 


92.3 Bi) Bey) 21.0 2.9 4.5 6.0 
Bye 217.2 28©=6 33. 1. 67.6 ~-18.6 26.5 10.7 
Be Cenlee zo. F52 -3.2 1 —.5 

-21.0 67.6 52 178.0 -37.0 -.9 83.1 all) 
2.9 -18.6 -3.2 -37.0 -27.6 -7.7. 10.0 
-4.5 DS it -.9 -7.7 48.9 4.6 
6.0  -10.7 -.5 83.1 1) .6 4.6 95.7 

FIGURE 13 


Figure 14 gives the numerical values of [¢] and the values of the 
seven natural frequoncics, Oe in radians per second, [6! 1s arranged 
in Figure 14 so that each eigenvector Ooh is listed below its related 


natural frequency Wee 


EIGENVALUES woos x103 


Zope ICO 2 TIS & 247.1 SpAlicge, 146.7. VsisAe 


BIGENVECTORS IN COLUMNS [9] 


17.7 Bee eon an 3 6 418.9 
=k BAG LY 2.8 iG 9 418.7 
mee 7s 9.5 SB Se hg 
-15.4 2 -1.2 Des y] 2.6 -627.3 | xl0o-' 
1535 44 -4.4 12.7 ile =o.) eae 
é 4.2 5 ane, oon Spo ese. s 
lone 1 2.5 ae 4 8.0 627.3 | 
FIGURE 14 








In the solution of the equations of motion by the technique shown 
in Appendix B, the time increment, AT, was 1 x noe o second. This value 
was chosen so as to be approximately one fourth the period of the mode 
having the highest frequency. If a value of AT was chosen close to one 
of the natural periods, the solution would become unstable and quickly 
diverge. 

Figures 15 through 32 show the results of different problems run on 
the model. The problems were designed to test each of the four configu- 
rations in which the model could exist. The cam profile was assumed to 
be harmonic with no dwells, and 1" peak to peak. There are four basic 
types of operating situations. 

i, eee eintained between cam and pushrod at all times without 


considering valve seat. In this situation the model could only exist 


2. Cam speed increased so that the pushrod no longer remains in 
ercact at all times with the cam. No valve seat. In this case the 
model will exist in Configuration I and II. 

3. Valve seat employed in the system. The cam rotated at a slow 
enough speed so that the pushrod would remain in contact with cam if the 
valve seat were not in the system. The model utilizes all four 
configurations. 

4, A combination of 2 and 3 above. The valve seat is employed and 
the cam is rum at a high enough speed for the pushrod to lose contact 
with the cam. Again all four configurations are employed in the model. 

The Figures 15 through 32 are graphs of the displacement and force 
curves at the contact coordinates. There are two graphs for each problem 


in Section A and B. The first graph will show the cam profile ara the 
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Be 











pushrod profile. The second graph will be the curve of F,, the force 


i! 
between cam and pushrod. Included on this graph for those problems 
employing a valve seat will be the curve of F., the force between valve 
and valve seat. The only deviation fron this will be the problem run 

at 20000 RPM without a valve seat. In this instance the pushrod profile 
was plotted on the same graph as the force profile. This was done to 
shay the comparison of dynamic responses of force and displacement at 
one coordinate. For section C, Figures 33 through 41, the graphs are 
plotted to show the most interesting features in terms of force and/or 
displacement of the particular problem. It should be noted that the 
values plotted are done at intervals of one degree which is much larger 
than the integration interval, AT, This gives a good representation of 
the overall force and displacement responses with one exception - 
MeieemrorcS Values. impact forses ccocur for cnly cne imtocration 


interval, AT, and therefore do not always get plotted. 


A. WITHOUL VALVE SEAT 
ie S000 RPM — Figures 15 and 16 
Figure 15 shows the pusiwod-cam profile curve and Figure 16 the 
curve representing the force between the cam and pushrod. As can he 


seen, the value of F, remains negative. As has been stated previously, 


u 


unless the value of Fy changes signs from negative to positive, the 
pushrod-cam profile curves will be superimposed on each other. 
2. 9524 RPM - Figures 17 and 18 
In Figure 18 the force reaches the zero value indicating sepa- 
Faction of the pushred and cam: Again it is noted that the value of F 


@ary gO nO higher than zero for the cam cannot pull on the pushred. Trere- 


meme, the flat part of the force profile indicates the length C= Ube tit 
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cam and pushrod are not in contact. Tne cam and pushrod are not in 
contact for approximately 50 degrees. When the two come back in contact, 
the force oscillates and is damped out by the internal damping of the 
system. Figure 17 shows the cam and pushrod profiles to be superimposed, 
but in actuality they should be separated for the same period of time 
igo) the Fy curve of Figure 18 is equal to zero. The pushrod is lifted 
such a short distance fram the cam surface that it is not distinguishable 
in Figure 17. This separation of cam and pushrod will be shown more 
effectively in Figure 19 for 11000 RPM. 
pm O00 RPM — Figures 19 and 20 

To campare with 9524 RPM and show a general trend of mechanism 
response the cam speed is increased to 11000 RPM. The pushrod gets 
‘thrown from the cam for a larger fraction of the cycle, about 180 decrees, 
emer Oure 20, and is thimewn furthor on Figure 19. Tf shouted Slee 


be noticed that the two curves are not superimposed on each other. This 


(D 


also happens in the 9524 RPM situation, but the separation was not great 
enough to be distinguished. The third feature of this problem is that 
upon coming back in contact with the cam, the pushrod bounces on the 
cam and again loses contact for a period before it finally comes in 
contact to remain. This is indicated by the separation of the pushrod-cam 
profile curves and also by the flat Omelet Tells F) curve. 

The response of the system to ever increasing cam speed can now be 
predicted. As the speed increases the pushrod will be thrown off the 
cam surface and the mechanism will remain in Configuration II, free of 
external forces, for longer periods of time. Also, as the speed 


increases, more and more bounces of increasing magnitude take place. Ihr 


addition, the impact force of collision between the two will in. 
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9524 RPM NO VALVE SEAT 


FIGURE le. 
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4. 20000 RPM - Figures 21 and 22 
In this case it can be seen that once the pushrod loses contact 
with the cam it never does come back in contact to remain. What takes 
place 1S a camtinuous bouncing as the pushisod and cam came in contact 
and the pushrod is throw off again. The period of time that the two 
do remain in contact is the time of take up of system inertia before the 


pushrod is thrown off again. 


B. WITH VALVE SEAT INCLUDED 

Next a valve seat will be included in the system. Upon impact, a 
and do are solved by the Soaniite te explained in Chapter IV, and the force 
between the valve and valve seat, Fo, is evaluated by Equation 27. a 


takes one time step, AT, to reach a value of 0. a4 takes one time step 


to reach the value (q. = Gg. ) /AT and enother time step to reach the 
/ 


C5 Cy 


value of zero. This means that Bo at time t, will usually be higher than 


the value at t. or t) - Again this is to be expected at the moment of 
impact. 
1. 1000 RPM ~ Figures 23 and 24 

As sham in Figures 17 and 18, until the cam reaches about 
9500 RPM, the pushrod will remain in contact with the cam if the model 
does not include a valve seat. Therefore, at 1000 RPM with a valve 
seat, any loss of contact between the cam and pushrod will be caused by 
the valve seat interference and not as a result of excessive cam speed. 

Figure 23 represents the response of the pushrod to the valve 
seat interference. Figure 24 gives the force at coordinate one, Fl, 
and the force at coordinate seven, Bo . On Figure 23 the cam and pushrod 


profiles are superimposed until the vaive comes in contact wit. oc) 
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seat. This contact is transferred through the mechanism and the pushrod 


is pulled fram the cam surface. In Figure 24 the upper curve is of F 


ab 
and the lower curve is of Fo. Notice that during the period of contact 
at the valve only, F., oscillates and Fy is equal to zero. ‘The reverse 


takes place when the cam comes back in contact with the pushrod. Fy 
oscillates and Fo goes to zero. The force curves contrasted against 
each other give the best representation of the particular point or 
points of contact of the mechanism. When either force curve equals zero 
Potsernaicates no contact at that particular coordinate. In addition, 
as will be shown for faster cam speeds, bounce of both the valve and 
pushrod takes place and is identifiable on the force curves. 

fee 2000 RPM - Figure 25 and 26 
This is a faster version of the previous example. ‘The force 
are of interest due to the oscillation upon impact. In this 


case, as evidenced by the additional flat areas on the Fy and FP, curves, 
there are bounces of the pushrod and valve after impact. Also of 
interest in Figure 25 is abe indication of the increased oscillation of 
the displacement of the pushrod after separation frem the cam. 
fee 2000 RPM — Figures 27 and 28 

Everything stated previously is even more in evidence at this 
Speed. Figure 28 indicates several bounces of the pushrod and valve 
after impact. The oscillation of q, on Figure 27, after the pushrod 


loses contact, is more pronounced. 
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pet 20> RPM = Figures 29 and 30 
The speed is above 9500 RPM so it can be expected that the valve 


seat and excess speed will have interacting effects on the model. Of 
interest 1S wnat happens first when the system is floating free of 
external forces - will the cam com@ in contact with the pushrod before 
the valve comes in contact with the valve seat? 

The pushrod first loses contact with the cam due to the action 
of the valve on the valve seat. The cam and pushrod come back in contact, 
the force Fy oscillates, and the pushrod bounces and is finally thrown 
free. At just before 460 deyrees, as shown by Figure 30, the valve comes 
in contact with the valve seat. Just after 480 degrees the pushrod also 
comes in contact with the cam and the system is subsequently thrown free 
‘again. In this case when the system is floating free of external forces 
the valve seat comes Lack in contact before tne pusnrcd. 

eee 20000 RPM - Figures 3] and 32 

The pushrod again first loses contact due to the valve seat 
interaction, comes back in contact with the cam, and is thrown free again. 
From then on the system resembles that shown on Figure 21 since the valve 
does not come in contact with the valve seat again. The pushrod continues 
to bounce on the cam at such a height that the valve does not reach the 


position of the valve seat. 


C. RELATED PROBLEMS OF INTEREST 
The previous examoles have shown the more routine results obtainable 
with the model. Figures 33 through 41 are graphs of sam of the additional 


Situations imposed on the model for examination. 
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1. Machining Error 

Figure 33 represents what happens due to a small machining error 
on the cam surface. The error was beteeen 268 degrees and 272 degrees. 
The value of cam displacement was assumed constant and the cam velocity 
and acceleration taken equal to zero. The best indication of the results: 
of the flat space is given by the oscillation of the force profile. On 
all graphs data points were plotted every degree. For this reason, this 
graph does not indicate the complete response. When the pushrod comes 
up to or goes off the flat space it hops off the cam surface for a small 
period of time, about .75 degrees, and then comes back to the cam surface. 
This should ke indicated on Figure 33 by the force going to zero. Since 
this occurs ketween successive data points it does not get plotted on 
the graph. The oscillation of the force at: coordinate one on Figure 33 
indicates probiens that could be associated with several machining errors 
on the same cam surface. The machining error here is 0.0003 inches. 

2. Undamped at the Lowest Natural Frequencies - Figure 34 

Figure 34 is the force between the cam and pushrod, Pos for the 
model without damping. This graph substantiates the value of the lowest 
natural frequency, 691.radians/sec. or 6602.5 RPM, of the undammeed model. 
This value from Figure 14 was obtained by solution of Equation 14 for 
the eigenvalue problem. At the resonance frequency of 6602.5 RPM it 
takes very little force above the static value to cause the mechanism 
to operate. As can be seen on Figure 34 the force oscillates about 
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Be olvyoyne Cam Profile = Figures 35 and 36 

Up to this point all the cam profiles have been harmonic, with 
no dwell, and represented by .5 sin (wt) inches with w equal to the 
value of cam speed in radians per second. Polydyne cam profiles are 
quite popular and this next example examines different responses of the 
mechanism to polydyne profiles. 

There are several references that are gcod on development of 
memyayne Cam profiles, (4), (10), (11), and (12). It is not the intent 
here to develop this theory, but just to outline how the profile was 
set up for solution on the model. 

Basically all that is needed for model input is q), = 6 , and 
qq: With this in mim a polynomial equation, 

om, * ©,° + C,0* + C,0° + C,0% + c.0° + C,0° +: C07 (29) 
was written with 0 equal to the angle of aoeeicn of the cam in degrees. 


Differentiating for oon qe and q gives, 


— A +2 3 ly 5 6 

ec, * 20.0 + 3C,0° + AC ,0 + 5C,0% + 6C 0 + 7C20° , (30) 
i aD Sas Eee 5 

i, = 2a 6C.,0 + 12C,0% + 20C,0° 4 30C,0°% 4 42C,0° , (31) 
eee _ - Zz x nm 3 > yy a) 
qd) = 6C., + 24C,0 + 60C.0* 4 120C,6 + 210C,0% . (22) 


The initial conditions were taken to be ay a and dy equal to 
zero at the beginning and end of the stroke. The value of q, at the 
beginning and end of the stroke would be specified. The one restriction 
placed on this situation was that the pushrod was not allowed to leave 
the surface of the cam at any time regardless of the value of the force 


at coordinate: one. 





Taking © equal to zero as the beginning of the stroke, Che oN 
and C3 are seen” to be zero due to the initial values of Gy ay and Sh 
Co equals the initial value of gq, and thus Equations 29 through 32 may 


be written in matrix form. 


4 5 6 fi 
qo 6) © Q Q C, 
: 3 4 5 6 
dy 40 50 60° 70 C, 
a (33) 
gq 1202 2003 300" 49° Ee 


ad 
OY 


a 240 6002. «12003 )~=3— 21094 c, 


The unknown coefficients, Cyr Ce, C., and C_ can now be solved 


G 7 
for by substituting =a and qq = Ay =a, = 0, the values at the 
F 


end of the stroke, into Equation 33. 


aly ~ <5 


0 


where dy |p indicates the value of q, at the end of the stroke. Using 
Equation 34 and the value of © at the end of the stroke in Equation 33, 


Cy Cc can be found. With Cyr Ce Cer and C.. known, the 


aay 7 7 


and C 
values of Gy qs and ay can now be solved for at any time within the 


stroke. A similar solution for a 345 polydyne profile can be made. 
However, only the initial and final values of Qy1 dys and q, nees ae 


defined. 
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Figure 35 shows the basic cam profile for both the 4567 profile 
and the 345 profile. Figure 35 also gives the response of the valve 
displacement q- to each of these profiles. Both cam profiles were taken 
to open in 90 degrees, to dwell for 90 degrees, to close in 90 degrees, 
and again to dwell for 90 degrees. This, in effect, simulates the 
inclusion of another cam at coordinate one so that the pushrod follows 
precisely the profile of the cam at all times, no matter what the value 
of F, might be. 

Figure 36 gives an example of the response of the systen when 
the amount of dwell is changed. A 4567 profile is simulated. The first 
run has a 90, 90, 90, 90 degree profile as previously described, but the 
second run takes 45 degrees to open, dwell for 45 degrees, close in 
225 degrees and dwells for 45 degrees. Again this figure gives both 
cam and valve profiles. 

4, 15000 RPM - Figures 37 to 41 

It is of interest to examine the relative deflections of indi- 
Vidual members of eS meithemtign fox the harmonic cam previously used, 
See Sin wt. 

Figures 37 to 39 represent the amount of deformation undergone 
by the model members during one particular operation. Figure 37 repre- 
sents the axial deformation of the valve stem. Figure 39 indicates the 
amount of bending undergone by the rocker arm. 


Figures 40 and 41 are the force and displacement graphs for 


15000 RPM in the same manner as shown previously. 
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re re 7. -_ 


On Figure 41 at about 850 degrees it is shown that both the cam 
and valve seat have come in contact with the mechanism. Careful exami- 
nation shows that the cam contacts prior to the valve. In the length 
of time it takes for the system to take up the inertia and be thrown off 
the cam again, the valve also comes in contact with the valve seat. It 
should be noted that 15000 RPM is between the 13205 RPM value of cam 
speed where it was noted in Figure 29 that the valve hits the valve seat 
before the pushred hits the cam and the 20000 RPM value of cam specd 
where the value never comes back in contact with the valve seat. As can 
be seen in Figure 41 for 15000 RPM scmetimes the valve does come in 
contact, i.e., 850 degrees and at other times it doe not, i.e., 500 


degrees and 1250 degrees. 
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VL. CONCLUSIONS 

An elastic link model has been developad and the equation of motion 
redefined to fit the external application of forces at two coordinates. 
This model has been used to study the gross effectsof action taken by 
the model as it is fitted into different loading situations. The solu- 
tion of the equation of motion was carried out by an integration 
technique. 

One of the problems associated with the model as now defined is that 
it does not take into account the propagation velocity of elastic waves 
in the structure. This gives some internal disturbances due to the fact 
that sane of the hich frequency mode shapes travel faster through this 
mathematical model than they would through an actual mcdel. This is best 
evidenced by tne initial increase in force between the cam and pushrod 
when the valve strikes the valve seat. A computer printed output of the 
action at coordinate one at the moment of impact at coordinate seven shovs 
that there is at first a reduction in Bis but that subsequently, before 
the pushrod is pulled away and the force goes to zero, it increases. An 
examination of Figures 23 to 32 indicates en initial increase in the 
negative value of the Ey curve before it is brought to zero with the 
separation of the cam and pushrod. This is due to the events occurring at 
one end being "felt" at the other end of the mechanism in one time ste. 
The high frequency mode shapes will reach peak value in the model prior 
to pushrod-cam separation. This occurrence is slight and happens in a short 
period of time, (depending on the speed, but it is measured in milliseconds) 
and the overall effect of the pushrod lifting from can surface is stiil 


not affected. 
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Another area of interest, which was not examined, was that of 
accounting for all parts of the model separating from each other. This 
would involve internal forces and would probably best be solved by 
working with the element, vice the system coordinates. 

The effect of clearance can be incorporated along the same lines as 
mentioned above. As of now the only clearences that can be incorporated 
are at coordinate one and seven. 

As the model now stands it gives a good indication of the effect of 
impact, machine error, cam speed, and valve placement on the action of 
the different members in the system. It could be argued that some of 
the forces abtained on impact would rapidly tear the mechanism apart. 
There is no counter argument to this except to state that the intention 
is to show the cause and inmediate effect in terms of force experienced 
by the system. it shouvuid be noted that ail runs are at a very high 
speed and therefore, suodject the model to very high forces. Interpre- 
tation of the large forces on impact, the oscillations of force, and 
the amount of bounce of one element on another, is part of the reason 
for setting up such a model. 

Many of the Gesign criteria for a good valve train pee een violated 
in order to demonstrate the potential of the model. The model is well 
adapted to examine features that should not occur in a good valve train 
and graphically denonstrates what will happen in terms of force and 
displacement responses if they do occur. 

The model is also extremely useful in examining a system that appears 
to be eating sufficiently well. Force oscillation in high speed 


mechanisms can be a problem which will not normally become evident until 


re) 





material fatigue appears. The mechanism does not have to run improperly 
memhave large force oscillation. Figure 33 for the machining error 
gives a good indication of this. 

Miers felt that the model in its) presente Gonniguration Can give an 
overall picture of system response. With relatively slight modification 
the model can be made to represent a specific mechanism and give 


mcividual response Characteristics. 
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APPENDIX A 


MATRIX MANIPULATION OF EQUATIONS OF MOTION FOR 
DIFFERENT SITUATIONS OF PROBLEM 


A. CONTACT AT CAM AND CAM FOLLOWER ONLY 


In the basic equations of motion 


im, aj} + fe, .1{4,} + [kj .}{a,} = (Fy) (35) 


wae |,2,...7) and (j = 1,2,...7). 


Fy 1s unknown and B es Fo are known and equal to zero. qy7 Gyr 
and qs equal the value of displacement, velocity and acceleration of 
the cam. 


Equation 35 can be expanded aid partitioned to cake aavantagc cr 


the known values. 


Expanding and partitioning Equation 35 yields, 
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= : (36) 


Removing the top equation for Fy , Equation 36, the set of remaining 


equations can be written, 





al A ey, do 
q, + 
Tal ee dy 
( . | 
“21 SIP Oy, 12 
a : qy 46 : : : (37) 
ea | 7208 a7 “y 
a oe ae Io 0 
=P qt oF = 
“oa Np a a 2 


The full set of equations can be rewritten as Equation 38 and 39. 


(m,.aq. + C,.d. + K459,) (38) 


a 
H 
toms 
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a ea - im] {q,} 7 {e),}4, + [ed] {q,} 


(39) 
+ Oky Jay + Uy )iq.} = (0) 
een — 2,5,...,/) (kK = 2,3,86¢,7) + 
Rearranged, Equation 39 can be rewritten as, 
where 
{p;} = = ea qd, 7 {c,13 qd, 7 Cae qy (41) 


Same 2,5,.+.,/) (k = 2,3,..-;/). 
This set of equations can now be solved by the numerical integration 
technique outlined in Appendix B, to yield iq,} : (q,), (q,}, where 


(5 = 2,3,...,7). With these values and knowing on and q,, Fj can be 


B. NO CONTACT 
In this case all the forces are known to equal zero. The basic 


equation of motion is left as is, 


im; 1{q, 3 - le; Jia; 3 + ik; 5) fay) = (0) (42) 


Paeee (i = 1,2,...,/) ( = 1,2,...,-/). 
Equation 42 is solved by the technique of Appendix B to yield {qty 


{q,) and {q,} where i = 1,2,...7. 
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C. CONTACT AT THE CAM AND VALVE SEAT 
The equation (35) is partitioned one step farther than for part A, 
for in this case Fo is not known, but qa equals the value of the valve 


seat and qd amd qy equal zero. Again as in part A, F, - Fe are known 


Le 
mercoual zero, Fy is unknown, and dy qq and qy equal the values 
associated with the cam profile. Equations 43, 44, 45, and 46 are 


formed in the same manner as part A, i.e., expanding, partitioning, 


and rearranging Equation 35. 


i 
Py SE sh Fh + ea! ey 
e ne @ 
By SE Ss, * i (44) 
fe 7 + [c. ey eee ic et (45) 
[ 3p) (953 ae hes) | 5x! 14 5 


fiere ] = 2,3,...6 and k = 2,3,...6; 


and {P.} - ENE ie) - ee Kee ey = on ol = {ko}, -- {e5)}qy 
Again the matrix equations can be solved by the technique of 
Appendix B to find (3) a5}, and te 3s (ee oe Os Fy and F> can 


then be solved by Equations 43 ana 44. 
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D. CONTACT BETWEEN VALVE AND VALVE SEAT ONLY 
In this case Fo is not known, qq and q, are egual to zero and qq 
is egual to the valve position. Expanding, partitioning, and rearranging 


Equation 35 yields Equations 47, 48, and 49. 


zs 
li 


O54 pee he Te 4a) (47) 


[fms 


mice Raves 7 Seis. 48 
im. tig. } a fess Isla + kod {q.3 (P5} (48) 


mgere (yj = 1,2,..-,6) (kK = 1,2,...,6). 


and ue }= tm, 7}dy - {c.5385 - {kha (49) 


Appendix B is used to solve the matrix equations and find {aj}, (a3, 
feet, 
4) J 


i 


1,2,...,6). F. is then solved by Equation 47. 
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| APPENDIX B 


INTEGRATION ALGORITHM 


A. DEVELOPMENT 


This algorithm is developed to solve the equation of motion, 


IM]{q} + [Cl{q} + [Kl{q) = {F) (51) 
by using a step by step time integration technique. This involves a 
choice of a function to represent how the acceleration changes within 
the time step T. 
Assume that, a, (t) = £. (t) (52) 
where 1 = 1 to the degree of freedom. 
Ss ies Greer 


Integrate for velocity and displacement; 


1Q 
oe 
i 


f£ (tat + Cy, = g(t) +O), (53) 


. = fg, (t) + Ct +O. = h, (t) + Ct +C,. . 
(54) 


ce) 
© 
| 


Express q, (t), q, (t) and q, (t) in terms of q. (t) evaluated at 
adjacent time step end points. 
sur Aa Aue Aiiane Sear Sein nln 
Pas P=2 Tale t ees 
The undetermined coefficients of £, (t) are expressed in terms of 
acceleration values at time step end points. The results are substituted 
into the equation of motion. This yields an algorithm which relates 


values of acceleration at adjacent time steps. 
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Assume a linear acceleration function within each time step. This 


yields an adaptation of Newnark's method of integration to matrix algebra. 








pape: 
z oe 
7 a 
di 
T ate 
n ntl 
mage n 
. ae “1H (qs = q.) 
Acceleration = qs (t) = eer ae E (55) 
(att ae D 
Velocity =< Be eer ee ee 
Y=q,(t) =q, +q.t + aR 5 (56) 
ess 
| heh “yn 2 (ap - ay) te 
Displacement = q, (t) = q, t+ qtt+q, 5+ mai aa (7) 
Evaluate qs {t) and qs (1) cites tee tee 
ar on en A eee 
qi es oe Ce 
ml ny en ap a Gt (wT) * 2 oot (ar) (59) 
4500 cat 4,73 4; 6 7 
n on ome ah 
= ¢. 60 
let a, =q. +4, 5 (60) 
n ae) n (aT)? 
= 61 
and yp. q, + q, 40 + gq, 3 (61) 
Now ihe and qs can be rewritten: 
ntl 7 ae 4 wnt ee (62) 





se ‘1 ie ae 
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2 
eel ntl “ny (AT) 
a2 13 Se 


qd. = Wy 


i i aah 6 (63) 


qd: nt+1 


Substitute the above results for an and qs sbaudey Weil SicsMein 5)1)- 


igh} + cleo) + Ftch ia) + oxi") + Eng egrtty = gph, 


(64) 


. As agent 
Rearranged, ts + SIC) + ney SSS = yee ee 


(65) 
. = Ns —nt+1 
Rewritten, (eskicn oe = {Pp ae (66) 
ideal 
where {p ~} = {F} evaluated at time ntl. 
my = im) + Sicj + Sy ‘a7 
2 6 

ad Beh ieee n n 

{Po} = fp = eer a eiearm (68) 


B. ALGORITHM 


ie cep 1 
The initial conditions can be determined from Equation 51 


evaluated at time equal to zero, 
my ige} + (c)iq2s + [KI{(q?} = (F9s 


Any two of {q°}, (4°) or (q?} may be specified and the third solved for 


by the above equation. 
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p. step Z 


Substitute the initial condition values into equations 60 and 61, 


Z 
Wet. pall n (AT) ze! 
yh tg a ae a 


eee) 


Set up fo fron Equation 67, 


Ra Z 
i = OP + Frc) + Fy 


Note that IM], tom and IK’ matrices may depend on motion values and 
this jay require modification during analysis. For linear providis, 
@ n why y 
1.e., this model, [M'], [C ], and [K’] remain constant and [M] is only 
computed once. 

4, Step 4 


Solve for {P} from Equation 68. 


_ ea n, _ n 
(Pi } = {py} ~ (Cltoe} ~ Ik tyy) 


De otep 5 


S08 
From Equation 66 solve for (q, Z : 
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GO ocep 6 


Evaluate Gt and {qh} from Equations 62 and 63. 


ae} = tar tig 4 
ont] nt+1 Picea (Gkal, : 
Moe ei ian ee 


ieee wtep / 


Evaluate {a} ana {y°'*; by Equations 60 and 62. 


ntl, love ma ronan uk 
nt]. (AT)? 


ml n+] ieee 
Cy} = tg} + {gh yar + {gy} 


Oo) 
9 


Step 6 
Let mn equal ntl. Go to step 3 and reiterate for 


“n+l 


(q" ntl, 


}, {g""} and tq 
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